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ABSTRACT 

We use the "Millennium Simulation" to study the mass function of accreted sub-halos 
during merger events in the dark halo assembly history. Our study includes three kinds 
of sub-halo mergers: (1) mergers that happen to the main progenitor of dark halos; (2) 
mergers that happen on the entire merging history tree of dark halos; and (3) mergers 
that leave identifiable sub-halos in present-day dark halos. We estimate the unevolved 
sub-halo mass functions (USMFs), for which sub-halo masses are measured at the 
times of their accretion. For sub-halos that merge into the main branch of a present- 
day dark halo, their USMF can be well described by a universal functional form, in 
excellent agreement with previous results. The same conclusion can also be reached for 
the USMF of all progenitors that have merged to become sub-halos during the entire 
halo merging history. In both cases, the USMFs are also independent of the redshift 
of host halos. Due to tidal disruption, only a small fraction of the accreted halos 
survive as sub-halos identifiable in the present-day dark halos. In cluster-sized halos, 
about 30% of the survived sub-halos are sub-subhalos, and this fraction decreases with 
decreasing halo mass. For given halo and sub-halo masses, the accretion time has very 
broad distribution, but the survived sub-halos are all accreted quite recently. 
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1 INTRODUCTION 

In the past two decades, the Cold Dark Matter (CDM) 
model of structure formation has been widely adopted and 
serves well as a framework for modeling the galaxy forma- 
tion. The understanding of CDM halo formation is a key 
factor to understand the formation process and properties 
of large-scale structure as well as the galaxies that form and 
evolve within the dark halos. There are multiple important 
properties regarding dark halo formation that have been 
intensively studied, including halo mass function (e.g . 
Bond et all Il99ll : lLacev fc Cole I Il993l: Istieth fc Torman 
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center of and merge along with dark halos, these events are 
therefore highly correlated with galaxy evolution. In this sce- 
nario, mergers play a transitional role in converting central 
galaxies into satellite galaxies in the post-merger halos. They 
may also trigger the evolution of various galactic properties, 
such as the morphology, luminosity , color and spacial distri- 
bution (|Barnes fc Hernquist I [1999; iNaab fc Burkert I [20031; 
Hopkins et al.ll2 006: Mailer ct al. 2006: De Lucia fc Blaizotl 



20071 : iMcIntosh et al.|[200a ) 



momentum property (e.g. | Barnes fc Efstathiou I Il987l : 
ICole fc Lacev I Il996l: IBullock et alj l2001a|). mass accretion 
histo ry (e.g. lZhao et al.ll2003l : IWechsler et al.|[20o3: iLi et all 
20071). and clustering pr opert y 
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Besides the properties mentioned above, halo-halo 
mergers and the subsequent evolution of resultant sub-halos 
have been of great interests rece ntly, in both analytical 
models and A-body simulations (Sheth! 120031 : ICao et ail 
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Although, from a phenomenological point of view, 
some observational statistics based on current sub- 
halo mass model matches well with the observations 
l|Mandelbaum et al.ll2006l : lKim et~aT1l2008l ). our understand- 
ing of galaxy formation still needs improvements due to the 
insufficient modeling of various physical processes, such as 
cooling, feedback, and merging history. It is also ambiguous 
how exactly post-merger galaxies are linked with pre-merger 
dark halos, because once a merger happens, the subsequent 
tidal forces and dynamical friction will cause the sub-halo, 
formerly a host halo, to loose mass and possibly become 
completely destroyed. This process gives rise to several pos- 
sible fates of the stellar compo nents of the gala xies that 
merge along with the sub-halos l|Yang et al.l 2009). 
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Despite the details of how a satellite galaxy evolves in 
a denser environment, it is always important to quantify 
the mass function of the associated sub-halos at the time 
of merging, for several reasons. First, previous studies have 
suggested that mass is a key fact or of various properties of 
dark halos, such as density pro file l|Navarro. Frenk fc White 1 
119971 ) , su b-halo populat ion ( Gao et al.l 12004 ). clustering 
property (|Mo fc White I Il99rj )~ Secondly, and more im- 
portantly perhaps, different approaches such as the halo 
occupation distribution (HOD) or similar models (e.g. 
Berlind fc Weinberg 1 12003 ; IZheng et al]|2005l ; iTinker et all 
20051 : IWang et all l2006l)~ and conditional luminosity func- 



tion (CLF) model l|Yang et al.ll2003l ; Ivan den Bosch il2007l ). 
which base their galaxy statistics on host halo mass, have 
resulted in reliable descriptions of the distribution of galax- 
ies. Therefore, to better understand the link between sub- 
halo mergers and post-merger galaxies, some important 
issues need to be addressed. For example, based on ex- 
tended Pr ess-Schechter formalism and direct iV-body sim- 
ulations, Ivan den Bosch et al l (|2005f) and iGiocoli et al.l 
(|2008ah found that the unevolved sub-halo mass function 
(USMF) of the progenitors that merged into halo main 
branch follows a universal form. Their findings are useful 
because the results can be linked to the number of central 
galaxies that may have turned into satellite galaxies through 
direct merger into a final halo. However, this information is 
insufficient to account for all incidences of mergers during 
the entire galaxy assembly history, because the hierarchical 
nature of CDM model suggests that sub-halos were indepen- 
dent host halos before the time of accretion, and it is likely 
they inherit the generic sub-halo population by the time 
when they became sub-halos. Thus, to investigate the pos- 
sible effects on the statistics of galaxy properties from the 
angle of sub-halos, one needs to further clarify two ques- 
tions. First, is the USMF really generic (i.e., does it depend 
on other quantities such as redshift than halo mass)? Sec- 
ond, what may be the difference if one takes into account 
the inherited sub-halo statistics of a sub-halo itself? Today, 
high resolution TV-body simulations provide a direct way to 
measure the merger statistics of dark halos to relatively high 
redshift with a good mass resolution. 

In this paper, we take advantage of a large iV-body 
simulation and its distinguished sub-halo statistics to answer 
the questions mentioned above. This paper is organized as 
follows. In section[2]we give a brief overview of the simulation 
and the algorithm used to construct the halo merging tree. In 
section[3]we describe in detail how to use the merging tree to 
identify halo-halo mergers during the halo accretion history, 
and further derive the unevolved mass function of the sub- 
halos characterized in several different ways. In section|4]we 
study the accretion time of sub-halos and mass function of 
sub-halos accreted at given redshift. Lastly in section [5] we 
summarize our results. 



2 THE SIMULATION 

In this paper we use the "Millenniu m Simulation" (MS ) 
carried out by the Virgo Consortium (|Springel et alj|2005h . 
We have used the same simulation data in an earlier pa- 
per to stu dy the age de pendence of dark halo spacial dis- 
tribution (|Li et alj fioOS). This simulation follows the evo- 



lution of 2160 dark matter particles in a cubic box of 
500 /i~ 1 Mpc on a side, with a mass resolution of approxi- 
mately 8.6 x 10 8 /i _1 Mq per particle. The simulation adopts 
a fiat ACDM model with f2 M = Sldm + ilb = 0.205 + 0.045 = 
0.25, where f2dm and fib stand for the current densities 
of dark matter and baryons respectively; the linear r.m.s. 
density fluctuation in a sphere of an 8 /i _1 Mpc radius, as, 
equals 0.9; and the dimensionless Hubble expansion param- 
eter h = 0.73. There are 63 snapshot outputs between z — 
and z = 80, with a roughly even placement in ln(l + ,z) space. 
The standard Friends-Of-Friends (FOF) algorithm with a 
linking length parameter b\ — 0.2 is used to identify FOF 
dark halos. Only FOF halos with more than 20 particles 
are resolved. Based on the FOF catalogue, each FOF halo 
is then assigned a corresponding "virial mass", M^, so that 
the average density contrast between the "virial halo" and 
cosmic critical density p c is approximately 200. The value of 
Mh is slightly smaller (about 5 — 10% on average) than the 
total mass of the corresponding FOF halo, and in general Mh 
accounts for the mass in the central part. This definition is 
therefore less affected by the so-called linking-bridge prob- 
lem that can cause uncertainties in the halo mass. In what 
follows, we always use "virial mass" , Mu, as the halo mass. In 
order to ensure robustness and completeness of our sub-halo 
analysis, we only use sub-halos with masses above a mass 
limit Mum = 2 x 10 10 h~ 1 M (s . This mass li mit is slightly 
higher than the re-simulated halos used by IGiocoli et al.l 
(2008a), but the simulation volume of the MS allows us to 
use many more halos to gain better statistics. 

The halo merging trees in the MS are constructed 
on the basis of sub-halos. In each FOF group, self-bound 
sub-str uctures (sub-halos) a re further identified using SUB- 
FIND i|Springel et al.ll200ll ). with the largest "sub-halo" be- 
ing the "main" halo. A sub-halo 1 at redshift zi is considered 
a progenitor of another sub-halo 2 at zi (zi > Z2) if the ma- 
jority of its most bound particles are in sub-halo 2. In the 
literature, merging history trees based directly on FOF ha- 
los are widely employed to study the mass accretion history 
of dark halos. The "sub- halo" -based linking algorithm used 
in the MS, however, has special advantages over the FOF 
merging tree in the study of the evolution of sub-halos. By 
definition, this algorithm e nables a more clear-cu t history 
tracer of sub-halos (see e.g.. iFakhouri fc Ma II2008I ). 



3 UNEVOLVED SUB-HALO MASS 
FUNCTIONS 

Mergers are important events during the lifetime of a galaxy. 
If there exists a one-to-one correspondence between a cen- 
tral galaxy and a host halo, then all galaxies are initially 
central galaxies at some high redshift in the hierarchical 
scenario of structure formation. Subsequent halo mergers 
play an crucial role in galaxy evolution, in the sense that 
a central galaxy will be formally transferred into a satellite 
galaxy and perhaps evolve passively afterwards without a 
significant amount of star formation. As mentioned before, 
galaxy properties may be highly correlated with their host 
halo mass. Understanding the mass function of the progen- 
itors of the sub-halos at the time of accretion is a key step 
in understanding the formation and evolution of satellite 
galaxies. In what follows we will refer the mass function of 
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sub-halos at accretion as the the unevolved sub-halo mass 
function (USMF), which reflects the fact that the sub-halos 
at the times of accretion have not yet been processed by 
dynamical effects, such as tidal stripping. In the rest of this 
section, we will discuss the USMF of sub-halos in the fol- 
lowing three categories: 



(i) In Sub-section 13. II we focus on sub-halos on the main 
branch of the merging tree, i.e., progenitors that directly 
merge with the main pr ogenitors of d a rk hal os. The same 
case has been studied bv lGiocoli et al. | l|2008al ). 

(ii) In Sub- section 13.21 we include all sub-halos that have 
merged into the entire merging tree of a dark halo. 

(iii) Finally in Sub-section l3.3l we focus on sub-halos that 
are directly identifiable in the present-day halos, the so- 
called "survived sub-halos". 



3.1 Main branch sub-halos 

To construct the halo main branch, we start with the final 
halo at a given redshift Zh (in this paper, Zh — unless other- 
wise mentioned) , and trace its most massive progenitor (the 
main progenitor) in the adjacent snapshot at higher redshift. 
We then repeat this procedure for the main progenitor till 
the progenitor mass is too small to be resolved. During this 
procedure, we also search the indices of all other progenitors 
that have directly merged into the main progenitor. If a pro- 
genitor was an independent halo before merger, we register 
its mass as well as the redshift at which it was accreted (the 
redshift information will be used later to study the mass 
function of sub-halos at given accretion time). This method 
eliminates cases where progenitors were already sub-halos 
of other more massive progenitors at the time of merging. 
With the information collected in this way, we are able to 
construct the USMF of main branch sub-halos. The results 
are plotted in Fig. [T] for host halos of different masses (as 
indicated) . 

We adopt the same functional form proposed by 
iGiocoli et alj (l2008al ) to fit the simulation results. Given a 
final halo (host halo) mass Mh and a sub-halo mass at ac- 
cretion, Maco the USMF, F, is written as 



'M a 



M h 



dN 



dln(M acc /A/ h 

' M, ^ 6 



{ mtJ exp - c Vw) 



(1) 



where N stands for the number of sub-halos that were ac- 
creted and a, b, c, d are fitting parameters. At the low-mass 
end (M acc /Mh — ► 0), this is a power-law, while at high-mass 
end (M a cc/Mh — > 1), the function decreases exponentially 
with (M a cc/Mh) . If the other parameters are fixed, a rep- 
resents the overall amplitude, b indicates the low-mass end 
power-index, c indicates the transitional point where the 
curve changes its shape, and al determines the steepness of 
the exponential decline. However, different combinations of 
parameters can result in F with similar shapes within the 
mass range probed here (log 10 [A/ a cc/Mh] £ [—4,0]). There- 
fore we do not intend to fit the result for each host-halo 
mass bin separately. Instead, we use all the mass bins to 
obtain an overall fit, which is shown as the solid line in 
the last panel. For comparison the best-fit USMF obtained 
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Figure 1. The USMF of main branch sub-halos. The upper two 
panels and the lower left panel show the USMF of 2h = ha- 
los with M h = 10 121 , 10 131 , 10 14 - 5 HMjj, respectively. Data 
points are the average over all halos with mass Mh, error bars 
represent the standard error of the average. For reference, in each 
pane l we also plot, with id entical dashed lines, the best-fit USMF 
from IGiocoli et al] l|2008ah . In the lower right panel, we summa- 
rize all the data from previous panels, and plot equation (fT} (in 
thick solid line) with an empirical set of parameters (as indicated 
in the panel) which provides a universal fit to all of our data. 



from IGiocoli et all (|2008al ) is shown in each panel of Fig. Q] 
as the dashed curve. 

Our res ults show an overall excellent agreement with 
the result bv lGiocoli et alj (|2008al ). The values of the fitting 
par ameters we obtain ar e very close to what were proposed 
by IGiocoli et alj l|2008ah . with only slight difference. For in- 
stance, we find the low-mass end power-law index b = —0.76, 
which is a slightly shallower than their —0.8. Our slope is 
chosen so as to reconcile the slightly higher "shoulder" found 
in the mass range log 10 (Af a cc/Mh) £ [—1-5, —0.5]. Note that 
for halos with Mh = 10 12 ' 1 Ii~ 1 Mq, we do not have data 
points that cover far enough into the power-law part. 

We also estimate the USMF for host halos identified at 
redshift zu = 1, and the result is shown in Fig. [2] Although 
the cosmic density field has evolved significantly during the 
time interval from Zh = 1 to 2h = 0, the USMF at Zh = 1 
has the same form as that for Zh = halos. All these suggest 
that the USMF of the main branch sub-halos has a universal 
form, independent of host halo mass and the redshift at 
which the host halo is identified. 



3.2 All sub-halos on the merging tree 

The merging history of a dark halo is in general quite com- 
plex. At lower redshifts, after a halo has assembled its main 
body, mergers may primarily happen on the main branch. 
However, at higher redshift, when a large fraction of the 
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Figure 2. The USMFs of main branch sub-halos for host halos 
at Zj, = 1. Different symbols represent the data points for host 
halos with different masses (as indicated), and the solid line is 
the universal fits we have obtained from Zh = halos. 



final halo mass was still part of smaller progenitors, merg- 
ers that take place on the sub-branches of the merging tree 
can no longer be neglected. In addition, the sub-halos that 
merge into the sub-branches may still present at the time 
when their host halos merge into larger halos. Although it is 
likely that most sub-halos that merge at high redshift may 
have already been dissolved by dynamical friction and tidal 
stripping by the time when the final halo assembles, the 
satellite galaxies that merged along with them may be more 
resistant to these dynamical effects. Therefore, it is interest- 
ing to investigate the statistical properties of these merging 
events. 

In order to quantify the USMF of all sub-halos in the 
entire merging tree of a halo, we start from the final host 
halo and trace back to all its progenitors that have ever 
merged as a sub-halo, regardless whether the merger takes 
place on the main branch or sub-branches. Once we found a 
merger between two independent halos, we register the mass 
of the sub-halo and the time of merger. 

Fig.[3]shows the USMF of all sub-halos in the halo merg- 
ing tree, in the same way as Fig. [T] for sub-halos in the main 
branch. Interestingly, equation {TJ still provides a good de- 
scription of the USMF in this case, although the fitting pa- 
rameters are different from those for the sub-halos in the 
main branch (see the solid line in the lower right panel and 
the values of the fitting parameters listed in the panel) . Com- 
paring the results here with those shown in Fig. [T] we see 
that the overall amplitude here is higher, due to the fact that 
sub-halos on sub-branches are also included. In addition, the 
increase in the amplitude is much larger for low-mass sub- 
halos than for massive ones, giving rise to a steeper power- 
law slope in the low-mass end - compare the data points 
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Figure 3. The USMF of all sub-halos that merged on the entire 
halo merging tree, plotted in the same way as in Fig. [T] Here, for 
reference purpose, the dashed lines in each panel represent our 
"universal" fit of the USMF of the main branch sub-halos (the 
same as the thick solid line in the lower right panel of Fig.[T}. Sim- 
ilar to Fig. [T] in the lower right panel we choose an empirical set 
of parameters (values as indicated in the panel) for equation JlJ 
and plot in thick solid line, so that it simultaneously fits all data 
points from the previous three panels. 



in each panel with the dashed curve that shows the fitting 
result of the USMF for sub-halos in the main branch. This is 
not difficult to understand. When we trace back in time to all 
branches on the merging tree, the number of sub-branches 
on the halo merging tree increases significantly with redshift 
due to bifurcation. Meanwhile, the average mass of progeni- 
tors drops dramatically because of mass conservation. Since 
our mass function is based on the unevolved merger progen- 
itors, more mergers of low-mass sub-halos are expected at 
higher redshift. 

In Fig. [4] we show the USMF of all sub-halos for host 
halos identified at z^ — 1. Clearly, this USMF shows remark- 
able agreement with that for Zh = host halos, indicating 
that the USMF of sub-halos in the entire merging tree also 
has a universal form. 

Let fivee and Fmb represent our universal fits to the 
USMFs of all sub-halos and the main branch sub-halos, re- 
spectively. Fig. [5] shows the ratio of these two functions as 
a function of Af acc /Mh. At the low- mass end Ftioc is about 
four times Fmb, while at the high-mass end they are nearly 
equal. The significant excessive rate of mergers at low-mass 
end seen in the ratio indicates the abundance of sub-halos 
that were accreted by the sub-branches of the merging tree 
(we will discuss in details later). These sub-halos may end 
up as the so-called s ub-subhalos when they finally settle in 
the main progenitor l|Yanget alj|2009h . 

As we have seen, the USMF of the main branch sub- 
halos is universal, independent of the mass and redshift of 
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Figure 4. The USMFs of all sub-halos for host halos at = 1. 
Same as in Fig. [2] different symbols represent the data points for 
host halos with different masses (as indicated), and the solid line 
is the universal fit we have obtained from Zh = halos. 




Loq 1n (M /M.) 

Figure 5. The solid line shows the USMF of all sub-halos divided 
by the USMF of main branch sub-halos, based on the two fitting 
results we have obtained. Dashed line is a reference line of y = 1. 



-I 1 r- 



all sub-halos 
sub°+sub 1 +sub 2 



sub (main branch) 



L°g, ( M ccc/ M h) 

Figure 6. Comparison between the USMF of main branch sub- 
halos (gray solid line, which is same as the solid line in the lower 
right panel of Fig. [TJ, all sub-halos (dark solid line, the same as 
the solid line in the lower right panel of Fig- EJ) , the sum of sub 
and sub 1 -halos (dashed line), as well as the sum of sub , sub 1 
and sub 2 -halos (short-dashed line) from the model prediction 
(equation [3]). 



host halos. This proposition has been adopted by some au- 
thors when modeling t he population of s atellite galaxies in 
dark matter halos (e.g. lYang et al.ll2009l ) . lYang et all (|2009h 
assumed that the USMF of main branch sub-halos is self- 
similar, and sub-halos can be divided into different "levels" . 
Since sub-halos can themselves be considered as host halos 
at the time of accretion, their sub-halos (referred to sub- 
sub-halos, or sub -halos) are also expected to obey the uni- 
versal USMF. Similarly, all levels of sub-halos (sub'-halos, 
i — 0, 1,2,3, ■ ■ ■, where superscript '0' stands for the main 
branch sub-halos) should have the same form of USMF. The 
summation of the USMFs at all levels should be equal to the 
USMF of sub-halos in the whole tree. To test this, we rewrite 
equation ([1} as 

diV 

n un ,o(M acc |M h ) = — — (2) 



— ( 



M acc 

"Mh" 



exp 



/ M acc y 



Since equation (J3JI is universal, it should apply to all sub 1 - 
halos (i — 0,1,2,3,- • ■). This allows us to calculate the 
conditional USMF of sub'-halos given the host halo mass 
M h , 



i(M acc ,i|M h ) 



Uun,o(M aC c,i|Macc) (M acc |M h )dM a , 



(3) 



Fig. [6] shows the comparison between the USMF of main 
branch sub-halos, all sub-halos, the sum of sub and sub 1 - 
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Figure 7. The "unevolved" mass function of sub 4 -halos, F' , . . 

sub A 

Different symbols represent different final host halo mass. The 
thick solid line is the universal form of the USMF of main branch 
sub-halos, -Fmb- 

halos, as well as the sum of sub , sub 1 and sub 2 -halos pre- 
dicted by equation ((3]). The best-fit parameters used in the 
calculation are indicated in the lower-right panels of Figs. [T] 
and [31 respectively. Clearly, as we include more levels of sub- 
halos, the summation of their USMF approaches asymptoti- 
cally to that of all sub-halos. It is also interesting that sub 1 - 
halos contribute the largest fraction of small sub-halos that 
are not included in the USMF of the main branch sub-halos. 
Note that in the mass range M acc /Mh G [—2, —0.3], the dif- 
ference between the dark solid line and the short-dashed line 
in Fig. [6] is more significant. It is unclear if this difference is 
real, or it is due to the limited statistics of the simulation 
data. 

3.3 Survived sub-halos 

In the two cases discussed above, the USMFs do not seem 
to depend on the final halo mass or redshift, and appear to 
be "universal" . However, once a sub-halo merges into a host 
halo, it will undergo a number of non-linear processes such 
as dynamical friction, which causes the sub-halo to merge 
into the center of the host, and tidal stripping, which causes 
it to lose mass or to be completely destroyed. Therefore, the 
number of survived sub-halos may be significantly lower than 
the sub-halo abundance described by the USMF. Note that 
there are two kinds of survived sub-halos: those that directly 
merged into the main branch, and those that were already 
a sub-halo of a larger progenitor when being accreted by 
the main progenitor. Throughout this paper, we refer to the 
former as sub A -halos, and the latter as sub s -halos, which 
are also known as sub-subhalos. 

After removing the sub s -halos from our survived sub- 
halo catalogue, we construct the "unevolved" mass func- 
tion of the sub A -halos. The quotation marks are used to 



indicate that a certain fraction of the main branch sub- 
halos have been completely destroyed, although the sub- 
halo mass used here is the mass at the time of accretion, 
M acc . Since the destroyed sub-halos are not included, we use 
F a ' ubA to distinguish this "unevolved" sub-halo mass func- 
tion from the USMF discussed previously. Fig. [7J shows the 
"unevolved" sub-halo mass function so defined for host ha- 
los with M h = 10 12 \ 10 13 - 1 , and lO 14 ' 5 /i _1 M , respectively. 
Apparently the shape of F' , A depends strongly on host halo 
mass. Unlike the USMFs discussed previously, for given host 
halo mass Mh, F' , A is not a monotonic decreasing function 
of sub-halo mass, but rather, its amplitude lowers when the 
sub-halo mass becomes very small. This is caused by the 
dynamical processes after the accretion of sub-halo. Note, 
however, that the value of M acc /Mh at which F' , A peaks 
depends on the mass limit Mn m adopted. There is a high 
probability that a sub-halo initially accreted with mass M acc 
slightly above Mn m to become smaller than Mu m during the 
post-accretion phase and thus to be marked as "destroyed" . 
In addition, smaller sub-halos are more difficult to survive, 
because on average they were accreted into their hosts ear- 
lier. 

Fig. [7J may be used to estimate the number fraction 
of sub-halos that survive the mass-loss process. For exam- 
ple, for halos with M h ~ 10 14 5 /i _1 M Q , about 62% of the 
accreted main branch sub-halos above the mass limit Mu m 
have been completely destroyed, this fraction increases to 
~ 78% and ~ 84% for M h = 10 13 ft _1 M and M h = 
10 12 /i _1 Mq halos, respectively. This trend may be under- 
stood since small systems start to accrete progenitors earlier, 
and so their main branch sub-halos are subject to mass loss 
and destruction for a longer time. The shape of F' . A for 
host halos at Zh = 1 is similar to that at Zh = 0. However, 
the similarity here is less meaningful, because the shape of F 1 
is highly affected by the non-linear effects d uring sub-halo 
merg ers, which is a very stochastic process (|Angulo et al.l 
l2008f) . 

The resu l t pres ented here is consistent with that of 
iGiocoli et al.l (|2008al . their Fig. 4), although their result is 
based on the evolved sub-halo mass function. They found 
that for small host halos, there are less sub-halos with the 
same fractional mass, M su b/Mh (where M su b is the current 
mass of survived sub-halos), than more massive host halos. 
In Fig. [7] we have showed that for smaller host halo mass, 
the amplitude of the "unevolved" sub A -halo mass function, 
_ F'„ kA , is also lower. In addition, according to lGiocoli et alj 
(2008a|, and reference therein), for evolved sub-halos mass 
function, the low-mass end is always higher than the high- 
mass end. Combining their results with our results of -F s ' ub A, 
we see that the majority of the smallest survived sub-halos 
are not the descendants of the smallest sub-halos initially ac- 
creted, but rather, the descendants of those that are several 
times more massive. Fig. [8] plots M acc , the mass at accre- 
tion, against the current mass, M su b of sub-halos, in 200 
host halos with M h « 10 136 h _1 M . Three different sym- 
bols denote different sub-halo accretion redshifts. It is clear 
that given M su b, sub-halos accreted earlier generally have 
higher A/ acc . At very low redshift (z £ [0,0.05]), sub-halos 
with a wide range of mass (log 10 [M acc /Afh] £ [—3.3, —0.5]) 
have been accreted by the main progenitor, and they have 
barely suffered from the mass loss so that M su b ~ M acc . 
However, for sub-halos that were accreted at high redshift 
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Figure 8. M acc against M su b for host halos of a given mass 
Mh. Points are from 200 randomly selected host halos with 
M h ft: 10 13 ' 6 /i _1 M s . Different symbols denote different redshift 
intervals during which the sub-halos enter the main progenitor of 
the host halo. Dashed lines indicate the mass limitation in our 
analysis, Mr lm . 
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Figure 9. The ratio of ^ s ' ubB t° ^ s ' u t,A ( see * ex * f° r details). 
Different symbols indicate different host halo mass. 



M h = 10 121 , 10 131 , and 10 14 5 /i _1 M Q , respectively. Clearly, 
a significant fraction of sub-structures were sub-subhalos. 



(z £ [1.2, 1.5]), their M acc are in general several times higher 
than M sub . 

Besides F' .a, we also construct the "unevolved" mass 

sub" ' 

function for the sub B -halo population, F' , B , in the same 
way as F' ubA . Fig. [9] shows the ratio of F' ubB to F 8 ' ubA , for 
host halos with the same masses as in Fig. [7] We would like 
to remind the reader, once again, that the sub-halo mass 
used here is measured at the time when they were last found 
as isolated halos. Given a sub-halo mass, the vertical axis in 
Fig. [9] is the ratio of the number of survived sub-halos ini- 
tially accreted by sub-branches to the number of survived 
sub-halos initially accreted by the main branch. In general, 
^sub B /^s'ub A * s hig ner f° r massive host halos. For a given 
host halo mass, though, this ratio is always low (~ 0.05) at 
the high-mass end (log 10 [M acc /Mh] > —0.7), because merg- 
ers involving sub-halos with mass comparable to that of the 
final host halo can only happen on the main progenitor at 
very late time. There also appears to be a generally increas- 
ing trend in this ratio as sub-halo mass decreases down to a 
certain point. This may be due to two reasons. First, some 
small sub-halos that merge to sub-branches of the merging 
tree may survive if the time scale for disruption is long. Sec- 
ond, as the redshift increases, the number of mergers that 
happen on sub-branches is not negligible. The increasing 
trend changes its sign when sub-halo mass becomes very 
small. The reason is that small sub-halos that are able to 
survive were most likely accreted in the recent past, when 
main branch already dominates the merger incidences. 

We can estimate the number fraction of sub s -halos 
among the whole survived sub-halo population, based on 
Fig.H This fraction is 9%, 17% and 28%, for host halo with 



4 ACCRETION TIME OF SUB-HALOS 

Although the USMFs give a quantitative description on the 
abundance of accreted sub-halos in the halo assembly his- 
tory, it does not include the time (redshift) when the accre- 
tion happens. In galaxy formation models, the epoch when 
central galaxies became satellites is crucial as the physical 
processes relevant to galaxy evolution after the merger are 
expected to be different. It is therefore important to incor- 
porate the sub-halo abundance at different redshift into our 
analysis. 



4.1 Sub-halo mass function at given accretion 
time 

4-1.1 Main branch sub-halos and all sub-halos 

We define the mass function of sub-halos at given accretion 
time (redshift) as follows, 



/(*) 



d£ _ _ 

dz dm(M 



dN(z) 



c/M h )dz' 



(4) 



where F is the USMF, M aC c and Mh stand for the mass 
of sub-halos at the time of accretion and the mass of final 
host halo, respectively. To obtain f(z), we choose a redshift 
interval Az around a given z, and only count the number, 
N(z), of sub-halos accreted during Az. Within the redshift 
range of interste, we found Az ~ 0.1 effectively eliminates 
the noise and result in a relatively smooth shape of /. 

Fig. [10] shows / of main branch sub-halos and all sub- 



halos for host halos with Mh 



10 1 



h 1 Mq. Interestingly, 
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Figure 10. The mass function at accretion, /, of main branch 
(upper panel) and all (lower panel) sub-halos, given accretion 
redshift z and host halo mass 

10 13.6 h-iM Q . Different symbols 
and lines represent the data points and their best fits according to 
equation JT) (with fixed 6 and d, see text for details), at different 
redshifts. 



in each case, / can still be described by equation |TJ) rea- 
sonably well. In addition, we found that the low-mass end 
power-index b of / are virtually independent of z, and is 
quite similar to the power-index we have obtained from the 
corresponding USMF. Since F = j fdz, it is expected that 
the integration of / over z reproduces the low-mass end 



power-index of F. The exponential shape of / (described by 
d) at the high mass end also shows no obvious dependence 
on 2. On the other hand, the amplitude of / and the tran- 
sitional point where / deviates from the power-law clearly 
depend on the redshift. By keeping b and d fixed at the val- 
ues obtained from the USMFs (b — —0.76, d = 3.2 for main 
branch sub-halos, and b = —0.91, d = 3.0 for all sub-halos), 
we fit / according to equation Jl]). Styled lines in Fig. llOl are 
the best-fits of / so obtained at the corresponding redshift. 

In Fig. 1111 we show the best-fit a and c against the red- 
shift z, for host halo with different masses. Panels on the 
left are best-fit a and c for main branch sub-halos, while 
panels on the right are best-fit a and c for all sub-halos. 
In general, a always decreases monotonically as z increases, 
which implies that more sub-halos are accreted at lower red- 
shift, especially for massive halos. Meanwhile, c shows pos- 
itive correlation with z, which means that, compared with 
small sub-halos, the number of massive sub-halos drops more 
q uickly as redshift incr eases. This disagrees with the result 
of iGiocoli etafl (|2008al ). Their Fig. 1 shows that the USMF 
of the main branch sub-halos accreted before the halo for- 
mation time Z{ is identical to the USMF of the sub-halos 
accreted after Zf, with proper adjustment in the amplitude 
a only. However, as we just mentioned, the number of mas- 
sive sub-halos drops more quickly at higher redshift, and 
therefore simply offsetting the USMF of high-redshift sub- 
halos along the vertical direction cannot reconcile the lack 
of massive sub-halos and reproduce the shape at high-mass 
end of the USMF. Our results suggest that the relative abun- 
dance of massive sub-halos becomes higher at low redshift, 
consistent with the hierarchical formation of dark halos in a 
CDM model. 



4-1-2 Survived sub-halos accreted by the main progenitor 

Given the time of merging, let us look at the mass function of 
sub-halos that survive as sub-structures in the final halo. We 
focus on sub A -halos, as the results for sub s -halos are similar. 
Based on the sub A -halo catalogue, we can register the time 
when they first became satellites of the main progenitor. We 
use / s ' ub A( 2 ) to indicate the same sub A -halo mass function 
at given accretion time defined in equation Q. 

Fig. Q2] shows / s ' ub A at different redshifts, for host halo 
mass M h = 10 13 ' 6 h~ 1 M Q . Note that the redshifts we used 
to plot / s ' u1> a(z) is, on average, lower than the redshifts used 
in Fig. 1101 because at higher redshift such as z > 1, f^ ub A 
becomes extremely small. Comparing Fig. [I2]with the upper 
panel of Fig. 1101 one can find both similarity and difference. 
At very low redshift (z = 0.07), / s ' ub A and / are similar, 
due to the fact that sub-halos accreted by the main progen- 
itor recently have a high survival rate. However, at higher 
redshift (z = 0.6), / s ' ub A becomes much lower than /, ow- 
ing to the dynamical effects that can effectively destroy the 
sub-halos accreted at early time. 

As we have shown in Fig. [7] the "unevolved" mass func- 
tion of sub A -halos, F' mhA , is not universal. Besides, the over- 
all amplitude of ^l' ub A) a ls° deviates substantially from the 
original USMF of main branch sub-halos, Fmb, especially 
at the low-mass end. The reason is clearly demonstrated in 
Fig. 1121 When redshift increases, / s ' ub A becomes increasingly 
lower, especially for small sub-halos. Since F' is the integra- 
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Figure 11. Best fit parameters a and c, given fixed b and d, of /, against redshift z. Panels on the left shows the result for main branch 
sub-halos, panels on the right shows the result for all sub-halos. Different symbols represent different host halo masses, as indicated in 
the figure. 



tion of /' over 2, it is therefore expected that F' . A would 
have the behavior shown in Fig. [7] 

4.2 Distribution of sub-halo accretion time 

In the previous sub-section, we have discussed the sub-halo 
mass function at accretion for given redshift. It clearly shows 
that the abundance of sub-halo accretion varies with red- 
shift. In general, more sub-halos were accreted at lower red- 
shift. It also seems that sub-halos with different masses may 
be accreted at different time. 

Given sub-halo mass fraction M acc /Mh and host halo 
mass Mh, Fig. [13] shows the number of sub-halo at the 
time of accretion as a function of redshift. Clearly, for fixed 
Macc/Mh, small systems start to accrete sub-halo earlier. 
For instance, dark halos with 10 12 ' 1 fe _1 MQ begin to acquire 
sub-halos with M acc = 0.03Mh at z — 5 ~ 6, while for halos 
with 10 14 ' 5 /i _1 Mq, this happens at z ~ 3. Compared with 
small sub-halos, large sub-halos enter the system fairly late. 



Nearly all sub-halos with mass M acc = 0.3Mh enter their 
host at redshift z < 1.5. 

For fixed host halo mass, large fraction of small sub- 
halos enter the system through sub-branches, especially at 
high redshift such as 2 > 1, while massive sub-halos (i.e., 
A/acc/Mh = 0.3) enter the systems only through the main 
branch, at relatively lower redshift. In addition, as discussed 
in Section [3] almost all survived sub-halos (sub^^-halos) 
were accreted at redshift 2 < 1, and more sub A,B -halos are 
likely to survive in massive systems. 

On average, sub-halo accretion of dark halos is de- 
termined by the initial CDM density power-spectrum and 
shows hierarchical signature. The sub-halo accretion for in- 
dividual dark halos, however, can be very stochastic. Let 
_Pi/2,A/ acc /A# h (z) denotes the probability distribution func- 
tion (PDF) of the redshift 2 by which the host halo has ac- 
quired 1/2 of the total number of the main branch sub-halos 
with fixed mass M acc /Mi,. Fig. 1141 shows P1/2, M acc /M h as a 
function of 2, for sub-halos with mass M acc = (2 — 5)%Mh. 
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Figure 13. Sub-halo mass function at the time of accretion against rcdshift z, given sub-halos mass M acc and final host halo mass 
Mjj. Each row represents one definition of sub-halo, and different columns represent different host halo masses, as indicated by the 
arrows. There are three lines in every panel. Think solid line is for sub-halo with M acc = 0.03M h , light solid line is for sub-halo with 
Ma CC = O.lMh and gray solid line is for sub-halo with M acc = 0.3M^. 



We choose M acc ~ 3% of Mh to ensure that the majority 
(> 75%) of our catalogued halos would have more than one 
sub-halo mergers with such sub-halo mass. Clearly, the red- 
shift covers a wide range: z G (0, 6), (0, 5), (0, 3) for halos 
with M h = 10 12 ' 1 , 10 13 ' 1 , and 10 14 5 /i _1 M , respectively. 
This indicates that, even for the same sub-halo mass and 
host halo mass, sub-halo merger is a highly stochastic pro- 
cess. 



5 SUMMARY 

Halo-halo merger is the basis of galaxy merger. The time 
of merger and the sub-halo mass at the time of merger are 
two important halo properties relevant for modeling galaxy 
formation. In this paper, we study the mass function and 
other properties of sub-halo mergers during the dark halo 
assembly history. We studied three kinds of sub-halos: main 
branch sub-halos, all sub-halos, and sub-halos that survived 
the dynamical disruption after merger. We also studied the 
redshift dependence and evolution of sub-halo mass func- 
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Figure 12. The sub -halo mass function at accretion, 



V 



given accretion redshift z and host halo mass 



10 13 ' 6 h 1 M.q. Different symbols connected with styled lines rep- 
resent results at different redshifts. 



mass (|Giocoli et al.ll2008al ). In addition, we found that this 
function is also independent of the redshift of the host halo. 

(ii) The unevolved mass function of all sub-halos that 
have been accreted during the entire halo assembly history 
is also a universal function that shows no host-halo-mass or 
redshift dependence. 

(iii) There are roughly the same or double number of sub- 
halos, with mass 1% or 0.1% of the final host halo mass, that 
were accreted by progenitors other than the main progenitor. 
The amount is significant considering the central galaxies 
that merge along with such sub-halos may be more resistant 
to dynamical disruption that destroy the sub-halos. 

(iv) The mass function of survived sub-halos at the time 
of merging is not universal, due to the fact that large frac- 
tion of sub-halos that merged at early time are destroyed by 
dynamical friction and tidal stripping. The fraction of sub- 
subhalos can account for up to 30% of the whole survived 
sub-halo catalogue in cluster-sized dark halos, and decreases 
with host halo mass. 

(v) In general, more sub-halos are accreted at lower red- 
shift. However, for given host halo and sub-halo mass, the 
accretion time has very broad distribution. Survived sub- 
halos are accreted late and therefore represent a very spe- 
cial subset of the total sub-halo population accreted into 
host halos. 




Figure 14. P1/2, M acc /M h ( 2 )> given main branch sub-halo mass 
-Mace = (2 — 5)%Mij. Different lines represent different host halo 



tion, as well as the distribution of the redshift at which a 
sub-halo is accreted. Our main findings can be summarized 
as follows: 

(i) We confirmed the previous result that the average 
unevolved mass function of main branch sub-halos fol- 
lows a universal functional form, regardless of host halo 
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